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THE CALCULATION OF EXPECTATIONS FOR CLASSES OF 
DIFFUSION PROCESSES BY LIE SYMMETRY METHODS 

By Mark Craddock and Kelly A. Lennox 

University of Technology, Sydney 

This paper uses Lie symmetry methods to calculate certain expec- 
tations for a large class of Ito diffusions. We show that if the problem 
has sufficient symmetry, then the problem of computing functionals 

of the form Ex{e ^'^^ g{Xs)ds^ reduced to evaluating a sin- 

gle integral of known functions. Given a drift / we determine the 
functions g for which the corresponding functional can be calculated 
by symmetry. Conversely, given g, we can determine precisely those 
drifts / for which the transition density and the functional may be 
computed by symmetry. Many examples are presented to illustrate 
the method. 

1. Introduction. Suppose that a stochastic process X = {Xt : t > 0} is a 
solution of tlie Ito SDE 

(1.1) dXt = a{Xt,t)dt + h{Xt,t)dWt, Xq = x. 

Here W = {Wt :t > 0} is a standard Wiener process. We would like to com- 
pute the expectations 

(1.2) 6,9(^i) = £;(e-^^'-/o ^(^^^-^^iXo = x). 

These types of functionals have numerous applications. For example, if A = 0, 
g{x) = X, then the expectation gives the price of a zero coupon bond in 
financial mathematics. If g{Xt) = fih{Xt) and h{Xt) does not depend on fi, 
then the expectation is the two dimensional Laplace transform of the joint 
density of {Xt, Jq h{Xs) ds). 

Most cases must be handled numerically, but recently considerable at- 
tention has been devoted to the problem of finding classes of diffusions for 
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which the transition density and expectations such as (1.2) may be deter- 
mined exactly. The papers [4, 5, 7, 8, 14] provide some details. 

By the Feynman-Kac formula, if the solution of the Cauchy problem 

ut = ax^<Uxx + f{x)ux - g{x)u, X > 0, 

(1.3) 

^(x, 0) = V3(x), rr G = [0, oo), 
is unique, then it is given by 

(1.4) ^x(x,t)=i5;,(e"/oS(^^)'^V(Xi)), 

where Xt satisfies the SDE dXt = f{Xt) dt + ^2aX^ dWt. We will use Lie 
symmetry methods to solve (1.3) by finding an explicit fundamental solution 
of the PDE as a function of the drift. We have a number of results which do 
this. The first is Theorem 3.1 below. 

The Cauchy problem (1.3) does not usually have a unique solution unless 
we impose certain extra conditions (see, e.g.. Part 4 of [18]). In fact, our 
methods actually yield more than one fundamental solution for the operator 
in (1.3). See Section 3. The selection of the correct fundamental solution for 
the purposes of computing the expectation must therefore be treated with 
care. 

The plan is as follows. In Section 2, for a given function (7, we show 
that the drift / must satisfy one of a family of Riccati equations if the 
PDE is to have nontrivial Lie symmetries. We also show how to explicitly 
solve the Riccati equations. In Section 3, we derive classes of generalized 
Laplace transforms of our fundamental solutions. Our methods require a 
time independent solution of the PDE and we show how to obtain these. 
Sections 4 and 5 contain examples and general formulae for the fundamental 
solutions. We also use group invariant solution methods to compute the 
precise form of the fundamental solutions in cases where the transforms are 
difficult to invert. 

In the final section, given a drift /, we determine those functions g for 
which the Cauchy problem (1.3) can be solved by symmetry and the expec- 
tation S,x^g{Xt) can computed. 

Each method for determining exact solutions has its own unique strengths 
and weaknesses. Specific techniques work well for particular class of prob- 
lems, but may fail when we move outside that class. For example, methods 
which use eigenfunction expansions and stochastic analysis have been ex- 
ploited successfully for many years. The literature is enormous, but see [14] 
for some recent work. However useful descriptions of the necessary eigenfunc- 
tions are only available for limited classes of PDEs and typically the eigen- 
function expansion of a solution cannot be explicitly summed. Eigenfunction 
methods also become more difficult when we move to higher dimensions. 
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Lie symmetry methods have natural advantages. They are easy to use, 
cover a broad range of cases and one can develop very simple tests to de- 
cide whether or not a PDE is amenable to symmetry analysis. The major 
limitation is that many important PDEs do not have any useful symmetries. 

However, for the rich class of PDEs amenable to Lie group analysis, the 
technique is extremely effective. We obtain explicit closed form expressions 
for the fundamental solutions in terms of the drift and do not need any 
changes of variables or measure. A major advantage is that the method can 
be extended to higher dimensions. An understanding of Lie techniques in 
one dimension is essential if we wish to apply them in higher dimensional 
problems. 

2. Calculating the Lie symmetries. A symmetry of a differential equa- 
tion is a transformation which maps solutions to solutions. In the 1880s 
Lie developed a technique for systematically determining all groups of point 
symmetries for systems of differential equations. -"^ Olver's book [15], is a fine 
modern account of Lie's theory of symmetry groups as are the texts by 
Bluman and Kumei [2], Hydon [12] and Ovsiannikov [16]. For a survey of 
some recent work on Lie symmetries, see [6]. 

In [7] Craddock and Lennox proved that if a PDE of the form ut = 
ax^Uxx + f{x)ux — fix^u has nontrivial symmetries, then there is always 
a symmetry arising from an action of 5'L(2,M) that is a classical integral 
transform of a fundamental solution. We generalize the results of Craddock 
and Lennox and apply them to problems in stochastic calculus. 

We will work with PDEs of the form (1.3). We observe that every PDE 
of the form (1.3) possesses time translation symmetries as well as scaling 
symmetries in u. That is, u{x,t + e) is a solution if u{x,t) is and so is 
cu{x,t) for any constant c. We term these trivial symmetries. We require 
that the PDE have nontrivial symmetries. 

Proposition 2.1. 7/7/2, the PDE (1.3) has a nontrivial Lie symme- 
try group if and only if f is a solution of one of the following families of 
drift equations. 



(2.2) 



(2.1) 



Lf = Ax'^-^ + B, 
Lf = Ax^-^ + 



8(2-7) 



3 



a 



.2 



where 



(2.3) Lf = ax 



.7 



( 



x^ f{x) 
M2-7) 



) 




x^-y{x) y 

2a(2-7)y 



+ g{x) + 



xg'{x) 



2-7 



^There also exist so called generalized symmetries, but we will not consider them. 
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7/7 = 2 then the PDE has a nontrivial Lie group of symmetries if and only 
if 

(2.4) Uf = A, 

(2.5) Uf = Alnx + B. 

With v{x) = ^^^^J'^^ we have 

TTf in \ , /(^) n \ ZM , xg'{x)\nx 
Uf = —v {^) + ^v ^^)"^ + 2 + 

Proof. The proof of this result is textbook, so we just summarize it for 
the 7 7^ 2 case. By Lie's method, we look for infinitesimal symmetries of the 
form 

(2.6) V = ^(x, t, u)dx + t(x, t, u)dt + (^(x, t, u)du, 

where we employ the convention = dx etc. We seek to determine con- 
ditions on and (p which guarantee that v generates a symmetry of 
the PDE. Standard arguments show that r can only depend on t and 
can only depend on x and t. Further, cj) must be linear in u. Lie's The- 
orem says that v generates a local group of symmetries if and only if 
pr^v[nt — {ax'^Uxx + f{x)ux — g{x)u)\ = whenever u is a solution of (1.3). 
Here pr^v is the second prolongation of v. (See Chapter 2 of [15] for the 
explicit prolongation formula for a vector field.) If 7 7^ 2, this leads to the 
conditions 

(2.7) i{x,t) = ^ ^ Tt + x'^/'^p, (l){x,t,u) = a{x,t)u + (3{x,t), 

2-7 



(2i 



= 2a(2-7)^ "" " W^t ^ " 

^1-7 

-f{x)Tt^ri. 



2<t(2-7)- 

Here r, r\ are functions of t only and satisfy 

_2;2-7 X^"'^/^ -(I-7) 

(2.9) - ^^^Pu + m = ^^ru - Lfn + Kfp, 

with Kf = ax^iix-r/'-^ - ^/(x))" + /(x)(|x^/2-i - ^/(x)). The 
function P{x,t) is an arbitrary solution of (1.3). The proof is completed 
when we compare the corresponding powers of x and recognize that r will 
only be a nontrivial function of t if / is a solution of one of the given 
equations. The 7 = 2 case is similar. □ 
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2.1. Solving the drift equations. If we set h{x) = f{x) then the ODE 
Lf = Ax'^~'^ + B \s equivalent to the Riccati equation 

(2.10) axh' -ah + + 2ax'^-^ g{x) = aAx^''^^ + loBx^'^ + C. 

Similarly for the second of this pair of drift equations. These Riccati equa- 
tions can be linearized by the change of variables h = 2axy' /y. 

We point out that it is natural to consider two different cases with equa- 
tion (2.10), namely A = 0,^4 7^ 0, because the Lie symmetries of the PDE 
take a quite different form if ^ = 0. 

Finding solutions of the drift equations in the 7 = 2 case is simplified 
by the substitution H{hix) = ^ — cilnx. Under which, for example, 
U f = A becomes 

a^H' -aH + ^H^ + 2ae9{0 = (2^ + f)^' + ^- 
Thus we have precisely the same form of Riccati equation. 

2.2. The 7 = 1 case. If 7 = 1 the drift equations are 

(2.11) axf -af + y^ + 2axg{x)=Ax + B, 

(2.12) axf -af + \f + 2axg{x) = \Ax^ + Bx + C, 

(2.13) axf -af + \f + 2axg{x) = ^Ax'^ + ^Bx^^^ + Cx- fa^ 

Here A,B,C are arbitrary constants. The factors of ^ and | multiplying 
A and B are a notational convenience. We will concentrate in this paper on 
the first two Riccati equations, which we treat separately. The third case is 
harder and will be discussed elsewhere. 

3. Applications of the symmetries. Symmetries map solutions to solu- 
tions. If the symmetry group is sufficiently rich, then one can construct 
extremely complex solutions of a given PDE from trivial solutions. We are 
able to obtain an integral transform of the fundamental solution by applying 
a symmetry to a suitable stationary solution. We begin with the A = case. 

Theorem 3.1. Suppose that h{x) = x^~^ f{x) is a solution of the Riccati 
equation 

(3.1) axh' -ah + \h^ + 2ax^-^ g{x) = 2a Ax^'^ + B. 
Then the PDE 

(3.2) ut = ax^Uxx + f{x)ux- g{x)u, 2; > 0, 
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has a symmetry of the form 



(1 + 4et)(i-7)/{2-7) " ^ i a{2- 7)2(1 + 4et) 
"^"p{^(^( (l + 4etW^) )-^("^ 



X u 



t 



(l + 4et)2/(2-7)' 1 + 4^7' 



where F'{x) = f{x)/x'^ and u is a solution of the PDF. That is. Us is a 
solution of (3.2) whenever u is. If u{x,t) = uo{x) with uq an analytic, sta- 
tionary solution then there is a fundamental solution p{t,x,y) of (3.2) such 
that 



00 

(3.3) / e-^y"'uo{y)p{t,x,y)dy = Ux{x,t). 

Jo 

Here Ux{x,t) = C/(i/4)a(2-7)2A- 

Equation (3.3) is a generalized Laplace transform. We may recover the 
fundamental solution by inverting the transform and dividing by uo{y). Be- 
fore proceeding to the proof we need a technical result. 

Proposition 3.2. The solution U\{x,t) in Theorem 3.1 is the Laplace 
transform of a distribution. 

Proof. This follows from the observation that Ux{x,t) can be written 
as a product of X'^ for some value and an analytic function G(l/A). Any 
function which is analytic in 1/A is a Laplace transform. Further X'^ is the 
Laplace transform of a distribution. The product of two Laplace transforms 
is a Laplace transform. Hence Ux{x, t) is a Laplace transform. See the table in 
[19], page 348 for the inverse Laplace transform of X'^ for different values of i' 
and Chapter 10 for general results on when a distribution can be represented 
as an inverse Laplace transform. □ 

Now we prove Theorem 3.1. 

Proof of Theorem 3.1. Lie's method shows that (3.2) has an in- 
finitesimal symmetry of the form 

/O.N 8xt ^ / 4x2-^ Ax^-^tf(x) ^ . . 2\ o 

3.4 V = + At^dt - — + '\' +(3t + AAt^ )udu, 

2-7 Vo"(2-7)^ cr(2-7) / 
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where (5 = '^^2--y' ■ Exponentiating this symmetry and applying it to a solu- 
tion u{x,t) yields C/e. The idea is that a solution U\{x,t) of the PDE should 
be obtained as 

/■oo 

Ux{x,t)= Ux{y,0)p{t,x,y)dy. 





The symmetry solution has the property that Ux{x,0) = e~^^^ ' uq{x). Thus 
we should have 

- ,,2-7 







''uo{y)p{t,x,y)dy = Ux{x,t). 



But this is the generalized Laplace transform of u^p. We thus need to show 
that Ux is a generalized Laplace transform of some distribution UQp. Since 

y uo{y)p{t,x,y)dy 

^l/(2-7)-l dz 



oo 



2-7 ' 

we must show that Ux is the Laplace transform of some distribution UQp. 
This follows from Proposition 3.2. 

To see that p is a fundamental solution of the PDE, observe that if we 
integrate a test function f{X) with sufficiently rapid decay against Ux then 
the function u{x,t) = Ux{x,t)(p{X) dX is a solution of (3.2). We also have 

u{x,0)= Uxix,0)^p{X)dX= uo{x)e-^'' '' ip{X) dX = uoix)^x), 



where <^ is the generalized Laplace transform of Next observe that by 
Fubini's Theorem 



oo roo - 

- ,2 — , 



UQ{y)^{y)p{t,x,y)dy= / / uo{y)(p{X)p{t,x,y)e ^ dXdy 

Jo Jo 

oo roo 2_ 

/ uo{y)ip{X)p{t,x,y)e~^y dydX 
Jo 

oo 

ip{X)Ux{x^ t) dx = u{x, t). 



We know that u{x, 0) = uq{x)^{x). Thus integrating initial data uq^ against 
p solves the Cauchy problem for (3.2), with this initial data. Hence p is a 
fundamental solution. □ 

Corollary 3.3. If g = in (1.3), then there is a fundamental solution 
p{t,x,y) with the property that 



-^y""'p{t,x,y)dy 
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4e{x^~^ + Aa{2 - j)H^) 



^^■^^ (1 + 4et){i-7)/(2-7) ^""P \ a{2 - 7)2(1 + Aet) 

"^^(^( (l + 4.tW7) )-^(") 
and p{t,x,y) dy = 1. Here e = \a{2 — ^)'^\. 

Proof. Since g{x) = we may take ^0(2;) = 1 in Theorem 3.1. Observe 
that C/o(3^) i) = 1- Thus p{t,x,y) dy = 1. □ 

For the 7 = 2 case we have a similar result. 

Theorem 3.4. Suppose that 

x'^ / f(x)lnx\" f(x) / / f(x)lnx\' a\ xlnxq'(x) , , 

Let ui^{x) he a stationary solution of 

(3.7) ut = crx'^Uxx + fix)ux - g{x)u 

analytic near zero. Then there is a fundamental solution of ( 3. 7) such that 

(3.8) r e-^'''^^^'^y^'p{t,x,y)u,{y)dy = Ue{x,t), 
Jo 

where 



(3.9) 



^.(x,t) = ^^exp(-^^(^"-)^ + ^--^/"- + (f-^ + -)--)^^n 
' ^ VI + 4et I a(l + 4et) J 



and F'{x) =e~^f{x). 

Proof. The proof is similar to the previous result. If the drift / satisfies 
(3.6), then the PDE has an infinitesimal symmetry of the form 

. . d . 2 ^ 
V = 4xtlnx— + 4r — 

(3.10) 

- + Inx^ f M _ iV + 2t + (4A + a)tA n 

\ a \ ax J J ou 

Exponentiating v and applying it to uq produces U^{x,t). This has the 
initial value ?7e(x,0) = uo(x)e~^^/'^^^''^^^ . The proof follows the same lines 
as previously. □ 
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3.1. Finding stationary solutions. For 7 7^ 2, to obtain a stationary so- 
lution uo{x) given h we need to solve ax'^Uxx + f{x)ux — g{x)u = 0. We 
divide by x'^~^ to rewrite this as axuxx + h{x)ux — x^~'^g{x)u = 0. If u{x) = 
'u(x)e-(^/(2<x))/ ^^g^ 2a^x'^u"{x) - {2a Ax"^-^ + B)u{x) = 0. Finally, 
the substitution z = x^~^ ^ u{x) = w{x'^~^) reduces this ODE to 

2a'^{2 - -ifz^w"{z) + 20-2(2 - 7)(1 - -i)zw'{z) - {2aAz + B)w{z) = 0. 

This equation has Bessel function solutions. Thus, given / and g we can 
always determine the integral transform (3.3) and hence express the funda- 
mental solution, at least up to a Laplace inversion integral. 

3.2. Finding different fundamental solutions. Different choices of uq will 
in general lead to different fundamental solutions. In fact we may some- 
times extract two fundamental solutions from a single choice of uq . We illus- 
trate the first situation with an example involving a squared Bessel process 
of dimension 3. The PDF ut = 2xuxx + 3tix has two stationary solutions, 
uo{x) = 1 and uo{x) = 1/y/x. The choice uo{x) = 1 leads to 

(3.11) / e-^ypit,x,y)dy= ,^ , exp^ 



(l + 2At)3/2" l + 2Ai 
which is easily inverted giving 

(3.12) p(t,x,y) = ^l=e-(^+J')/(2i)g.^^j^%^^_ 

This is the transition density for the squared Bessel process of dimension 3. 
If we use the second stationary solution we get 

From which we obtain a second fundamental solution 
(3.14) p2(t,x,y) = ^^e"(^+^)/(2*)coshf 

This is not a transition density. Whereas Jq° p{t,x,y) dy = 1 we have 



00 



(3.15) / P2{t, X, y) dy = J^e"^'^^'^ + erf 





The function (3.15) is a solution of ut = 2xuxx + 'iux which satisfies 
\mit^Qu{x,t) = 1 for all x > 0, but it is not continuous at the origin. All 
solutions generated by p2{t,x,y) have this feature. They solve the Cauchy 
problem for the PDF on any domain [e,oo), e > 0, but are not continuous 
at the origin. 
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Next consider the equation 

(3.16) ut = xuxx + ^3 - ^ Ux, x>0. 

V o + ax'^/ 

An application of Corollary 3.3 yields the existence of a fundamental solution 
p{t,x,y) satisfying 

(3.17) / e-^yp(t,x,y)dy = ^ ^exp<^-- 

^ ^ Jo ^ ' {b + ax^){l + \tf I 1 + AtJ 

We take a, b positive. To invert this Laplace transform requires the introduc- 
tion of distributions. (A complete discussion of the inversion of the Laplace 
transforms which require distribution theory is found in [9].) Inverting the 
Laplace transform produces the fundamental solution 

(3.18) 

+ t{b+ax^) ^^y^^b^^^^y^- 

Here 8'{y) is derivative of the Dirac delta function, I^, is the modified Bessel 
function of the first kind. It is easy to show that p{t,x,y) dy = 1. If 
q{t,x,y) =p{t,x,y) — ^^:;^8' {y)i then q is also a fundamental solution and 
/o°° 9(^5 ^) y)dy = 1 as well, since 6' acts by /q°° (j){y)6'{y) dy = — (/)'(0). Indeed 

we may also take r{t,x,y) = ^ b+ax'^ e~^^"'"^)/*/2( ^^^ ) and this is a funda- 
mental solution with 

(3.19) r r{t, x,y)dy = l- e""/* + = /(x, t). 

Jo t{b + ax^) 

Now Z(x, t) — > 1 as t — > for all x > 0. Thus (3.16) with the data n(x, 0) = 
1 does not have a unique solution. If we wish to interpret p{t,x,y) and 
q{t,x,y) as probabilities we are faced with a situation in which we have two 
densities arising from the same generator that define the same null set and 
probabilities, but because of the terms involving generalized function they 
do not define the same expectations. 

The drift function / has the property that —1 < /(x) < 3. This means that 
if we consider in a formal sense the SDE dXt = f{Xt) dt + \/2Xt dWt we are 
in a situation where Xt may become negative because of the drift and hence 
the diffusion Xf may then become complex valued. It is not obvious how to 
interpret this. The more natural SDE dXt = f{Xt) dt + y^^2\Xt\ dWt leads to 
the PDE Ut = 2\x\uxx + fix)ux, but our theory does not cover a PDE of this 
form, though one can also analyse PDEs of this type by symmetry. A full 
discussion of the problem of trying to associate a diffusion with a PDE of 
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the form (1.3) covered by our methods is beyond the scope of this paper. 
But it is clear that our techniques ahow us to construct very general classes 
of PDEs which have very complex behavior. 

For our purposes we will restrict attention to problems where the inter- 
pretation of the fundamental solutions as densities is clear. That is, we will 
look at generators which arise from well behaved SDEs with positive drifts 
and/or unique strong, real valued solutions. We will include examples of fun- 
damental solutions with delta function terms, but our examples will be more 
straightforward than the above. A full investigation of the phenomenon we 
have described here will be considered elsewhere. 

This leaves the question of which choice of uq we take in order to use the 
Feynman-Kac formula to compute expectations. The examples we give will 
show how this choice is made. 



Example 3.1. Consider the Bessel process X = {Xt : t > 0} where 

(3.20) dXt = -^dt + dWt, Xq = x. 

Xt 

Bessel processes are well studied and the transition density for this process 
is known. Here a = 1/2,7 = ^'^cl f{x) = ^. We will require a > ^ for sim- 
plicity. So h{x) = a and we easily see that ^ = in the Riccati equation. Set 
9{x) = We want a fundamental solution of 

1 a a 

ut = -Uxx + -Ux - -r^u, /i > 0. 

2 X 4j;^ 

— (/i/4) 

If we are to calculate Ex{e ° (l){Xt)) then we require a fundamental 

solution that reduces to the known transition density of the process at /x = 0. 
A stationary solution is given by solving 2x^n" -|- Aaxu' — fiu = 0. We take 

the solution ^0(3;) = x*^ where d = ^ — a + -|- (a — ^)^. Observe that as 

^ — > 0, uq{x) — > 1. Setting e = Xa = ^ in (3.3) gives 

(3.21) / e-^y yMt,x,y)dy= , ...^. expf 



[l + 2Xty "^V l + 2Ai 

with u = d + a + ^. This is a Laplace transform if we let y"^ = z. Inverting 
the transform yields 

(3.22) M<, ^, !/) = ? f exp ij^'" 



t\x) V 2t ) \t 

If = 0, this reduces to the well known transition density of a Bessel process. 
We can now write Exie * ^'^'^ ^-^o ^) = e~^yp{t,x,y) dy. 
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This gives the Laplace transform of the joint density of {Xt,jQ -y^)- We 
cannot evaluate this integral. However 

2\ (2i/-2a-l)/4 







-x2 

e 



T{a)iFi{a,u,x'^ /{2t + At'^X)) 

^ r(z.)(i + 2a)" ' 

where a= ^+^°+^'^ and iFi{a,b, z) is Kummer's confluent hypergeometric 
function. 



There is another solution x"^^ with (i2 = ^ — a— Y^ + (a — ^)^. We observe 

that as /i — > x'^^ — > x'^ where v = ^—a—\a—^\ which is nonconstant. So the 
fundamental solution coming from this stationary solution does not reduce 
to the necessary density. This situation is typical of these problems. 

Note. There are many Laplace transforms which can be computed by 
other methods. However our methods often yield different representations 
of these quantities than can be found in much of the literature. For a Bessel 
process of index ^ there is a well known alternative representation 

(3.23) 

^7-5 /-oo 7;P-iexp(-x2(t; + A)/(l + 2(u + A)t)) , 

dv. 



T{p) Jo (1 + 2{v + A)t)i+7 

where p = ^(7 — ^-^id 7^ = + /U^. This is proved using the absolute 

continuity law = (— )^e ^ ^ Px^'^lj^t- However we have never 

seen the integral (3.23) explicitly evaluated. Our representation allows for 
the explicit evaluation of the integral from tables. Moreover to obtain the 
joint density, we only have to invert a one dimensional Laplace transform 
rather than a two dimensional Laplace transform, since we already know the 
inverse Laplace transform in the A variable. 

Example 3.2. Now we turn to the Bessel process with drift. The SDE 

is 

dX, = (^^- + ^-(^ j dt + dW,, a > -I. 

These processes arise from the radial part of d-dimensional Brownian motion 
with drift and were studied by Pitman and Yor in [17]. Here cr = \ and 
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h = xf satisfies ^xh' -^h + lh^ = !^x^ + l{Aa^ - 1). So vl = ^5^. Using the 
stationary solution uq{x) = 1 we obtain 

(3.24) / e^^y p(t,x,y)dy= °^ '\ V exp ^ ^ — -]. 

^ ^ Jo ^' ^ {l + 2\t)Ia{hx) ^\ l + 2\t J 

We convert this to a Laplace transform by setting y'^ = z. Inverting the 
Laplace transform gives the transition density 

Since 7^^^ — > (f as 6 — > 0, this reduces to the transition density of a Bessel 

~XXt~fi f' '''' 

process as 6 — > 0. Now consider the problem of computing Ex{e " ). 

We require a stationary solution of the PDE 

1 , f a + 1/2 , bla+i{bx) \ fi 
(3.26) Ut = -u^x+{ \ ^ . ]ux 

2 \ X ia\bx) J x-^ 



We look for a solution of the form ^0(2;) = e ^^^^vix) where F' = /. We 
see that 

ix'^v'ix) + 4(1 - 4a^ - Ab^x^ - 8fi)v{x) = 0. 
Since F{x) = ln{^/x Ia{bx)), we obtain the stationary solution Uq(x) = 
I "(fax) — ■ ^ote that as fi ^ 0, tig ^ 1. Theorem 3.1 tells us that there 
is a fundamental solution q{t,x,y) of (3.26) such that 

roo 



(3.27) 



Uo{y)q{t,x,y)dy 



il+2tX)Ia{bx) 

Inversion gives the fundamental solution 

(3.28) q{t,x,y) = ^^e.J-^-^ - ^-t]l r^(^-^]. 

Clearly this reduces to the transition density at /x = 0. Obviously 

Ex{e ) = /q°° e yq{t, x, y) dy, but we have not been able to eval- 

uate this integral. However to find the joint density, we only need to invert 
the Laplace transform q(t,x,y) in ^u. This can be done numerically. 

The procedure we have is now clear. We begin with an SDE which pos- 
sesses a unique solution and compute a transition density for it. This can be 
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done by our methods. We choose the fundamental solution of (1.3) which 
reduces to this density when 5 — > 0. From this we may compute the desired 
expectations. 

4. Laplace transforms of joint densities. As is clear from Example 3.1, it 
is possible to use the techniques we have developed here to compute Laplace 
transforms of joint densities very easily. We will focus attention on the 7 = 1 
case and do an example of the 7 = case below. It should be clear that 
we may prove similar results for any 7. If 7 = l,g{x) = fix^, with n = ±1 
and n = ^, then we can find drift functions which do not depend on fi. It 
is possible to easily compute Laplace transforms of the joint densities of 
{Xt, Jq Xg), (X^j/q-^) and indeed we can also find the Laplace transform 

of the joint density of {Xf, X^ds, Jq -y-)- The four dimensional Laplace 
transforms of the law of {Xt , /q X, ds Jq^, Jq VXl ds) can be computed for 
certain classes of diffusions, but this is more complicated and will be treated 
elsewhere. 



Example 4.1. Let X = {Xt :t > 0} be a squared Bessel process, where 

dXt = ndt + 2^/YtdWt, Xo = x. 

_ r' ds 

To compute u{x,t) = E^{e ^Jo Xs h{Xt)), we require a fundamental solution 
for the PDE 

(4.1) Ut = 2xUxx + u. 

X 

This is the case when cr = 2,7 = l,g{x) = /i/x in Theorem 2.1. The stationary 
solutions are uq{x) = x"^ , d± = ^{2 — n i: \/{n — 2)^ + 8/^). We take d = dj^, 
as x'^- diverges near zero and does not tend to the constant solution as 
^ ^ 0. Applying the symmetry we have 

(4.2, ^%^„M..).-«*=^^^^i;^»p{^}. 

This Laplace transform may immediately be inverted to give 

N . N 1 /x\(i^"/2)/2 fy/^\ r (x + y)] 

(4.3) P(M.rf = ^(-) W-i(V)'='"'|- a }■ 

From this we obtain the two dimensional Laplace transform of the joint 
density of {Xt, Jq -^ds). 



-^yp{t,x,y)dy = e--/'' 



2t 



^T{a)iFi{a,(3,x/{2t + 4f\)) 
r(/3)(l + 2At)° '' 
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with a = d+f,/3 = 2d+§. 

We have not been able to invert this Laplace transform. It is interest- 
ing to compare this result when A = with the famous Hartman- Watson 
law for a squared Bessel process which gives the conditional expectation 

EAe-"-'''''l'*^\X,^v)= ''^^"\ where 1 - >. 
Example 4.2. Now we turn to the SDE 

(4.4) dXt= ^ .^ dt+V2XtdWt, Xo = x,a>0. 

1 + aJi.f 1 2 



For the PDE ut = xUxx + ~ have two stationary solutions 

ui{x) = '^2+aT ~ ^^2(2;) = — — ' neither of which is equal to 1 
at /i = 0. However, the linear combination 2u2 + aui does satisfy the desired 
condition, so we use the stationary solution uo{x) = x(^"V^+^^)/^( ^+^^^^^'^'^ ). 
Letting d-t = i(l it ^/T+'4Jl) we obtain 

, , 1 / ax^+ 2x'^- \ / Ax 



2 + ax V(l + At)2<i+ {l + Xt)^d-J 1 + At 

Inverting the Laplace transform yields the fundamental solution 

t 1 + ax\j y 2y"-+ + ay'*+ V * 



(4.6) + 



2 + ay 2x"'- 
(2 + ax) 2y'^+ + ay<^+ 

xe--(^^yy^J'lo[^y,{^-y)d^. 



Here = C-^[X'^^+^]. This inverse Laplace transform is a right sided 

distribution. See [19] for a table of these inverse Laplace transforms. We 
must be careful when taking the limit, however a rather laborious calculation 
shows that as ;U ^ this reduces to the fundamental solution 



(4.7) pit,x,y) 



{2 + ax)t 

which is the transition density for the process Xt- From (4.6) we may calcu- 
late Ex{e ' ^Jo Xs), The integral can be evaluated exactly, but we leave 



ft d. 


this to the reader. 



We may readily establish more general results. We present an example 
shortly. First however, we need a lemma. 
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Lemma 4.1. The PDE ut = crxu^x + f{x)ux — ^u, with f a solution of 
the Riccati equation crxf — af + ^/^ = Ax + B has a stationary solution 
Uq{x) which has the property that Uq{x) = 1. 

Proof. This is a direct calculation. Settine; f = 2ax— reduces the Ric- 

t> J y 

cati equation to 2ax'^y" — {Ax + B)y = which has general solution y = 
^/x{clIa{V'2'Ax) + C2Ka{V'^Ax)), wherc a = ^\/2B + o"^ and Ka is the usual 
modified Bessel function. Now a stationary solution of the PDE is given by 
solving axu" — f{x)u' — = 0. The substitution uq{x) = e~^^^^/^'^'^^v{x) 
leads to the ODE 2a'^x'^v" - {Ax + B + 2^ia)v = 0, where F'{x) = f{x)/x. 
Thus there is a stationary solution of the PDE 

(4 8) u^(x) - "^^-^"^^ + (^2K,{V2M) 

CiIa{V2A^) + C2K^{V2M)' 



where i/ = ^ ^25 + + 4//cr. Clearly = 1. □ 



Theorem 4.2. Suppose that f is a solution of the Riccati equation 

2- 



axf -af+y^ = Ax + B with A>0. Then the PDE 



Ut = axuxx + f{x)ux — —u, 2B + + 4/icj > 

X 



has fundamental solution 



p{t^x,y) = X^e--^('=)/2'^-('^+^)/(''*)-(^/(2'^)*) 



at u'^{y) 

(4.9) J^^dhmj^y 



at J \ a 



at J \ a 

where F'{x) = f{x)/x, u'^{x) is as above and v = ^ \/ 2B + cj^ + Afia satisfies 
\iy\ < 1. From which we may calculate E{e~^^'~^Jo xS) by 

(4.10) ^,(e-^^'-^/o^)= r e-^yp{t,x,y)dy. 

Jo 



Proof. Let f{x) = 2a xy' /y where y{x) = ^(ci/q,(^^^^) + C2I- 
(Since a is not an integer Ka.{z) = Ia{z).) To obtain a stationary solu- 
tion of the PDE we use as in the previous lemma the substitution uq{x) = 
e-^(^)/(2o")^(3;) in the equation axu" + f{x)u' — fi/xu = 0. Then v must sat- 
isfy 2ax'^v"{x) - {Ax + B + 2ai^i)v{x) = 0, which has solutions \/xI±u{^^^^^). 
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Set 

We then have Wo(y) = 1- Now 



\ <7 



(4.11) 



Uo{y)pit,x,y)dy 

^^-F{x)/{2a)-\{x+At'2 /2)/{l+\at) 

l + Xat 



We use the fact that 



lexp 



o 9 

m + n 
A 



2mn\ 



Id{2m^)Id{2n^), d>-l. 



Inverting the Laplace transform gives the fundamental solution. The expec- 
tation follows from the Feynman-Kac formula. □ 

For outside the range of the theorem, the inversion of the Laplace 
transform will in general involve right sided distributions. Techniques for 
inverting these transforms are discussed in [9]. 



5. The other Riccati equations. Suppose that the drift / is such that 



h- 



r.l-7 



/ satisfies a Riccati equation of the form 



axh' -ah + i/i^ + 2axg{x) = Ax^""^^ + Bx^'^ + C, A / 0, 

then we do not get Laplace type transforms of the fundamental solution, ex- 
cept in some special cases. Instead we get a more general integral transform 
of the fundamental solution, known as a Whittaker transform. A similar 
situation arises with the final Riccati equation, but we will consider that 
elsewhere as there are some technical issues which need to be carefully dis- 
cussed. 

The reason for the difference is that the Lie symmetries for these Ric- 
cati equations are different in form to the symmetries in the previous case. 
Consider the 7 = 1 case. Craddock and Lennox proved in [7] that if / is a 
solution of the second class of Riccati equations with A^O and uq (x) is a 
stationary solution of ut = crxuxx + f{x)ux — ^x^u, then so is 



Ue{x,t) 



-Bt/{2cr) 



Axe 



exp 



+ 



1 



(5.1) 



X (e 



VAt 



VAt 



F 



xe 



VAt 



oV~At 



F{x) 



■.B/(2crV~A) 

e) ^ 'uo 



xe' 



oVAt 
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In general this does not reduce to the necessary form for a Laplace trans- 
form when we set t = 0. Moreover, for certain stationary solutions, the action 
is trivial. However, there are special cases when it is a Laplace transform. 
We present an example. 

Example 5.1. Consider the PDE 

ax 
I + ax/2 

The drift is a solution of the Riccati equation axf — af + ^/^ + = 
^Ax^ + Bx + C, with A = 4:fi,B = C = 0. There are two stationary solu- 
tions, (x) = e^v^^(ax + 2)^^. The symmetry (5.1) leaves Uq unchanged. 
However for the other solution we have 

5.2 Ue{x,t) = -^expi-^^ . 

2 + ax \ e V — ej 

By our general argument, with A = we can find a fundamental solution 
p{t, X, y) such that 



2 + ay 



(5.3) 



2 + ax V A(e2v^* - 1) + 2^e'^^* J ' 
Inverting the Laplace transform yields the fundamental solution 
, 2 + ayey^y ( -2^{x + ye"^^^) ] 

P^'^'^y^ = 2T7^7v^'A Tvw^i 1 

fjlx h{2./JI^ / sinh(^t)) ^ ^^^^ 
y sinh(^t) 

Letting /i ^ gives the transition density for the process X with dXt = 
i4-aXt/2 + \/ dWt . Consequently we can compute the Laplace transform 
of the joint density of the processes {Xt,rQXsds). Formula 6.621.3 of [11] 
and the Feynman-Kac formula give 

rt ^-Jlix-2JTix/(e^VT^^-l) 



X 2 + 



2 + ax 

^^(g(/3V4a) (4/^ ^ al3/{2a^)) - 4//5) 



sinh(y7''t) 

Here a = A — '^'^J%i^^ and /3 = 2y^^cosech(y^t). 
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Although we have recovered the necessary fundamental solution by invert- 
ing a Laplace transform, for these classes of Riccati equations we in general 
have to deal with a transform involving Whittaker functions.^ 

Definition 5.1 (The Whittaker transform). The Whittaker transform 
of a suitable function (p is defined by 

(5.4) <I>(A)=/ {Xy)-''-'/'e-^y/'W,+y2A^y)cPiy)dy. 



An inversion theorem for this transform is known. For a suitable constant p 
we have, in the principal value sense 

<^(y) = _^ r(l + z.-fe) r+*°°(Ay)-fc-i/2eA,/2^^ .{Xy)HX) dX. 

The functions Wfc_|_i/2,iy and Mi,_Yj2,u are the Whittaker functions given 
in 13.1 of [1]. For a discussion of the transform, see [3], page 110. Special 
cases of the Whittaker transform include the Laplace transform and Hankel 
transforms. 

In [7] the following result is established for the 7 = 1, g{x) = [ix^ case. The 
proof of the following result is identical to the case in [7] . 

Theorem 5.2. Let f he a solution of the Riccati equation 
axf -af + + 2axg{x) = ^Ax"^ + Bx + C. 



Let U^^^^_^{x,t) be given by setting e — > — 1 in 

Ue{x,t) = e-(^t)/(2-)(eV^* - e)S/{WA)-/3/2^;3/2g-F(x)/(2<x) 



2(j(ev^*-e) J V ' 'cj(ev^* 
where ^'(a, /3, z) is the Tricomi confluent hypergeometric function and F'{x) = 
f{x)/x. Suppose that X^/^"^^U^i^^^-^_-y^{x,t) is the Whittaker transform 
of a function hf^{t,x,y). Then this function is given by 



where p^ is a fundamental solution of the PDE ut = crxuxx + f{x)ux — g{x)u. 



^ The transform is due to Meijer and is often called a Meijer transform. However 
there are several other integral transforms which are also called Meijer transforms in the 
literature and this nomenclature can be somewhat confusing. 
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Similar results can be proved for the general case when 7 7^ 1. The Whit- 
taker transform is analytic, thus one can in principle invert a Whittaker 
transform as a series expansion, since the Whittaker transform of is of 
the form CA"""*"^ where C is a constant depending on n and the parameters 
in the transform. Unfortunately, tables of Whittaker transforms tend to be 
very sparse, so in practice, these transforms are really only invertible in the 
special cases where they reduce to Laplace or Hankel transforms. Further, 
theorems guaranteeing that a given distribution is a Whittaker transform 
are available only for limited parameter ranges. At present therefore, these 
results are of mainly theoretical interest. 

There are some other issues that arise when we try to calculate expecta- 
tions of the form Ex{e~^'^^~^ Jo Xsds-^^ T^h.^ stationary solutions of the PDE 
ut = ax^Uxx + f{x)ux — f^xu do not necessarily have the property that a 
linear combination tends to the constant solution when /i ^ 0. See Example 
5.1. Conversely both of the linearly independent stationary solutions may 
converge to a constant when we take the limit. We thus need a different 
criterion for deciding which fundamental solution to use. 

An alternative approach to the construction of the necessary fundamental 
solutions is to use group invariant solutions. A group invariant solution is a 
solution which is left invariant under the action of a group transformation. 
Chapter three of Giver's book [15] contains a detailed discussion of group 
invariant solutions. See also the material in [6]. We will construct the gen- 
eral form of the fundamental solutions and determine the necessary one by 
examining the integr ability properties. 

The following result, which combines two theorems in Bluman and Kumei's 
text [2] is the key to our method. 

Theorem 5.3. Consider the nth order boundary value problem P{x, 
D'^u) = subject to the conditions Bj{x, u, n^"""*^^) =0 on the surface ujj{x) = 
0. A vector field v is admitted by the BVP if: 

(1) pr"v[P(a;,L>"n)] =0 when P{x,D°'u)=0. 

(2) v{ujj{x)) = when iOj{x) = 0. 

(3) pr^~^v[Bj{x,u,u^'^~^^)] = when Bj{x,u,u^'^~^'^) = on the surface 
ujj{x) = 0. 

Suppose that a BVP admits a vector field v. Then the solution of the BVP 
is a group invariant solution with respect to the symmetries generated by v. 

We prove several results as an illustration of the technique. We concen- 
trate on the case ^ > 0, leaving the case ^ < to the interested reader. 
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Theorem 5.4. Suppose that f is a solution of the Riccati equation 
axf - af + \ f'^ + 2a/ix^ = ^Ax^ + Bx + C,A> {). Then there is a fun- 
damental solution of the PDE 

Ut = aXUxx + f{x)Ux - jJLXU, x > 0, 

of the form 

P{t,x,y) 

, , ^/3^e-(n^)-F(2/))/(2<x) / Bt + \ 

(5.6) = 1= exp -= 

2cJsmh(\/Zt/2) V 2(t 2cjtanh(^/At/2) / 



where v = ^'^'^+'^'^ and we interpret I^y{z) to he Ky{z) if u is an integer. 

Proof. Lennox proved in [13] that if the drift / satisfies the given 
Riccati equation, then the PDE has a Lie algebra of symmetries spanned by 

VI = dt, V2 = xe-^'dx + ^^dt - ^ {yJx + f{x) + e-^'udu, 



„—VAt 1 / , R \ r- 

V3 = -xe-^'dx + -^dt - — VIx - f{x) + e-^'udu, 



lA 2a \ ' ' VA 

V4 = udu, v = P{x,t)du- 

Here /? is an arbitrary solution of the PDE. The symmetries correspond- 
ing to adding a solution /? to the original solution. Plainly these cannot give 
group invariant solutions. We will solve 

Ut = (TXUxx + f{x)Ux — fJ.XU, X > 0, 

(5.7) 

u{x, 0) = 6{x — y). 

We look for symmetries of the form v = X]fc=i CfcV^ which preserve the 
boundary conditions. We require the boundary x = to be preserved by the 
action of v. That is v(x) = when x = 0. We also require the boundary t = 
to be preserved, which means that v(t) = when t = 0. Finally u(x,0) = 
5{x — y) must be preserved, which requires v(m — 6) = when u(x,0) = 
5{x — y). These conditions are satisfied if and only if C2 = C3,ci = — -^C2 

and C4 = 2^^^^^^C2. This gives the form of the symmetry we need. 
We then compute the invariants for the action of 

/ 2 iAy + B) 

V = Ci ^Vi + V2 + V3 + 2 7^^V4 

V y/A ay/A 
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Invariants of the action are given by solving the equation v(r/) =0. They 
are readily found by the method of characteristics and we may take the 
invariants to be rj and v where 

(5.8) !]■ ^ 



(5.9) 



4sinh^(V^t/2) 
u = exp 



{Bt + F{x)-F{y)) ^/A{x + y) \ 
2a 2cjtanh(VAt/2)/ 



X V 



Asmh^{y/At/2)J' 



Now v{ri) must satisfy 2a'^T-fv" {ri) — (C + 2Ayri)v{r]) = if u is a solution of 
the PDE. Hence 

Since the solution of our initial value problem is a group invariant solution 
for the action generated by v the result follows from substitution of v and 
1] into the expression for u. □ 

Note. The symmetries used in the proof are actually symmetries for 
any PDE ut = axuxx + f{x)ux — g{x)u when axf — af + ^/^ + 2axg{x) = 
^Ax"^ + Bx + C. So we may use the same technique to construct the form 
of the fundamental solutions for any PDE of this type. We will present the 
case when g{x) = ^ below. 

5.1. The functions Ci{y),C2{y) and transition densities. The functions 
Ci, C2 will depend upon the initial and boundary conditions we may wish to 
impose on the PDE . The Bessel function I^y{y) is not integrable near zero 
for > 1, so if we require a fundamental solution which defines solutions 
of the Cauchy problem for constant initial data, then we will have C2 = 0. 
It is also this choice with Ci = 1 which is needed if we are to recover the 
transition density as we let /z ^ 0. In the general case, Ci, C2 can be found 
from the condition that u{x,t) = (p{y)p{t,x,y) dy solves the PDE with 
u{x,0) = (j){x). For most cases, one may simply take 4>{x) = l,(j){x) = x for 
the initial data, from which we may identify Ci and C2. 

Recovering a transition density needs care however. Though our expres- 
sion for the fundamental solutions are very general, it is possible to have 
fundamental solutions for the types of PDEs we have been considering which 
involve terms with delta functions. We have seen such a situation in Example 
5.1. 

In all save one of the examples we present below, we do indeed obtain the 
transition density by letting /i ^ in the fundamental solution we obtain 
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from Theorem 5.4. However in general, letting /i = will produce a fun- 
damental solution of ut = axuxx + fix)ux, but there is no guarantee that 
this fundamental solution will be the transition density for the diffusion. We 
need to be careful to incorporate the behavior at the boundary y = 0. 

One can usually produce a density from the fundamental solutions pro- 
vided by Theorem 5.4. Suppose that p^{t,x,y) is the fundamental solution 
we obtain from an application of Theorem 5.4. We then have 4>iy) x 
p^{t,x,y) dy — > cj){x) as t — > 0. Now introduce the function h^{t,x) = 
Pfj,(t,x,y) dy. If lim^^Q hf^{t,x) = 1, then will reduce to a density at 
H = 0, since an application of Lebesgue's dominated convergence theorem 
gives liuifj^^Q p ^{t, x,y) dy = 1. If linifj^^Q h^{t,x) / 1 then to produce a 
density, we let p{t,x,y) = lim^^oP^tCi, a;, y) and h{t,x,y) = limfj,^oh^{t,x). 
Then p{t, x, y) + {1 — h{t, x, y))6{y) is a density. 

In this case to calculate the expectations we want, we will have to incorpo- 
rate additional terms involving generalized functions to produce a fundamen- 
tal solution that reduces to the necessary density. Consider the symmetry 
solution (5.1), where we choose a stationary solution not left invariant by 
the symmetry. We make the following observation: lim^^o Ui{x,t) = 0. This 
follows from the fact that the stationary solutions of the PDE (5.7) are of 
the general form 

^X0(3;)=X^/2g-(F(x)-v^x)/(2a) 

(5.10) 

' ^ { ^ VAx\ ^ ( ^ ^^Ax^ 
ciiFi a,/3, +02^- a, /3, 



with a = ^(^ f7(l + 2C/cj^)), (i = ^J\^2C j a'K The estimates 

iFi(a, 6, z) = ^e^z'^-^(l + 0(\z\~'^)) and ^(a, 6, z) = z~''{\ + 0(\z\~'^)) for 
large \z\ with z > give the result. (See page 504 of [1].) Thus if = 
p^{t,x,y) + Ui{x,t)5{y), then 

oo ;>oo 

(p{y)q^{t,x,y)dy= / (l){y)p^{t,x,y) dy + <j){0)Ui{x,t) ^ (j){x) 
Jo 

as t — > 0. Hence q^ is also a fundamental solution. 

We apply this to the diffusion in Example 5.1. An application of Theorem 
5.4 leads to the fundamental solution 

, 2 + aye^y ( -2^{x + ye'^y^^)] 

^'^^''"'^^ = 2T^;7^^M ^^^1 1 

^ //Ix /i(2^//Ixy/ sinh(^t)) 
V y sinh(^t) 

If we add Ui{x,t)6{y) to this we obtain the necessary fundamental solution 
found in Example 5.1. We do not have a proof that this method works in all 
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cases, but similar procedures will always allow us to construct the necessary 
fundamental solutions so that we can compute the desired expectations. 

Corollary 5.5. For n>2 a fundamental solution of the PDE 

Ut = 2xUxx + — —u, X > 0, 

is 



p{t,x,y) 



2sinh(bt) \x 

(5.11) 



^ r 2tanh(6t) ; I sinh(6t) ) ' 

Proof. We simply need to check that / satisfies the Riccati equation 
and determine A,B,C. We choose Ci = 1,6*2 = 0. We further note that 
as b ^ 0, the fundamental solution converges to the transition density for 
an n dimensional squared Bessel process. One may find other fundamental 
solutions which do not reduce to the transition density. □ 

We now have a very easy proof of a well known result for squared Bessel 
processes, which is usually proved using martingale methods. 

Corollary 5.6. Let X = {Xt:t > 0} he a squared Bessel processes, 
where 

dXt = ndt + 2^/YtdWt, n > 2. 

Then 

i?,(e-^^'-(''/2)/o^^'^^) 

(5.12) 

_ exp(-(x6/2)(l + 2A6"icoth(6t))/(coth(6t) + 2A6"i)) 
~ (cosh(6t) + 2A6-1 sinh(6t))"/2 ' 

Proof. The expectation is given by 

(5.13) E,(e-^^'"('''/2)/o^.rf^) = r e-^yp{t,x,y)dy, 

Jo 

where 

p{t,x,y) 



n/4-1/2 



2sinh(6i) \x 

(5.14) 

f b{x + y) \ ( b^ \ 

^ """P r 2tanh(6t)j ^(""^^/^ l^hlhM j ' 
Evaluation of the integral is routine. See the tables in section 6.6 of [11]. □ 
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5.2. The zero coupon hand price in the CIR model. Using our results, it 
is easy to compute the well-known zero coupon bond price in the CIR model. 
We have A = h'^ + A^a, B = -ah, C = - aa. With Ci{y) = 1, (72(1/) = 
and F{x) = a\n.x — hx we recover the necessary fundamental solution. We 
leave this to the reader. Let us now consider an example where we have to 
include a delta function term. 

Example 5.2. We are interested in the process X where 



dXt = 2Xt tanh(Xi) dt + ^/2XtdWt, Xq = x. 

Using the symmetry (5.1), one may check that the transition density of this 
process is 

(5.15) ,(M,ri = 4^i^exp{-f±ii|(.&,(HVl|U,i„,,,,,) 

smh t cosh x { tanh t } \ \ y \ smh t ) 

Using our results we find a fundamental solution of 

(5.16) Ui = xua;a; + 2a;tanh(x)ua; — AiXU 
is 

1 coshy r k(x^y)\ \kx (2^fkxy 
^' ^) = sinh kt cosh X ^"Pr I V 

Here k = y/1 + /x. This does not reduce to the transition density as — > 0. 
However, uo{x) = e~^^~^f^^ sechx is a stationary solution of this PDE. If we 
apply the symmetry (5.1) to this solution, we find a new solution Ue{x,t) 
which has the property that Ui{x,t) = exp{- ^^^g^J . 

Now q{x, y, t) = p{x, y, t) + 6{y)Ui {x, t) is also a fundamental solution and 
this does reduce to the transition density. From which we may compute 

Ex(e~^^*"^^o^^'^") 

-^yq{x,y,t)dy 



00 

e~ 





„—xl tanh(fct) 

= \J\(x t) + - ^gfcxCsch(fct)/(cosh(H)+(-l+A)sinh(fct)) 

_|_ gfcxCsch(fet)/(cosh(fct)+(l+A)sinh(fct)) _ 2^ 

5.3. Calculating expectations of the form Ex{e * '^Jo-^s). In Section 3 
we calculated the expectation CA,/x/x(^t) for the case when 7=1 and crxf — 
o"/ + i/^ = Ax + B,A>0. We now show how to compute these expectations 
when / satisfies (2.11). 
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Theorem 5.7. Suppose that axf -af + \P = ^Ax^ + Bx + C, A>^. 
Then the PDE 

(5.18) ut = axuxx + f{x)ux - ^u, f^>0, 

has a fundamental solution of the form 
p{t,x,y) 



(5.19) 



Ae^Fiy)-F{=^))/{2a) 



/X r_^_^I(x+^| 

Vy I 2(T 2(jtanh(VAt/2) J 



2crsinh(\/At/2) My \ 2a 2atanh{y/At/2) 

''^^^^'^(.sin^!lt/2))+''^^^)'-^(.sin^!lt/2)))' 



in which F'{x) = f{x)/x and v = ^^^^ interpret I^y{z) to he 

Kv{z) if V is an integer. 

Proof. The proof is similar to the previous result. The invariants are 
the same, but the PDE in this case reduces to the ODE 

Aa'^rfv"{r]) - (2C + AAyr] + Afia)v(j]) = 0. □ 

The following result for Cox-Ingersoll-Ross process is an easy corollary. 

Corollary 5.8. The PDE ut = axuxx + (a - bx)ux - with /i > 0, 
a,b> has a fundamental solution 

P^^^^^y)- 2asmh{bt/2)\xJ 

^'■''^ ^^^K^("* + ("-^)-i^^niw2)) 

^ ''\asmh{bt/2)J' 
Here v = ^^J{a — cj)^ + 4/xcj. 

Corollary 5.9. Let 

k = ^, a = T^fl + coth^^^^ +A, P 



a & A , fbt\\ , ^ bJ^ 

— , a = — 1 + coth — +A, P= , 

2a 2a\ \2 J J ' 2cjsinh(6t/2 



and Ms^r{z) he the Whittaker functions of the first kind. For the CIR process 
dXt = (a - bXt) dt + ^/2aXt dWt we have 

^ ' r(i/ + i) 
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(5.21) X exp( 7^ ( at + X 



2a V tanh(6t/2) 



Proof. Observe that the fundamental solution in Corollary 5.8 reduces 
to the transition density of a CIR process as — > 0. So we have 

(5.22) ii;,(e-^^*-''/o^)= r e-^yp{t,x,y)dy. 

Jo 

The result follows from the fact that 

X^-/V«».,(.,^,...4±^^e-.-M_,,(^), 

which is formula 6.643.2 of [11]. □ 

5.4. Laplace transforms of joint densities for {Xt, Jq Xg ds, Jq By 
similar means to the above, we may find fundamental solutions of the PDE 

ut = CTXUxx + f{x)ux - (I + ^J^.x)u. 

Theorem 5.10. Suppose that f is a solution of the Riccati equation 
axf -af + + 2av + 2a^ix^ = \Ax^ + Bx + C,A>{). Then the PDE 

(5.23) Ut = axuxx + f{x)ux — i — \- fix]u, /.f>0, z^>0, 



has a fundamental solution of the form 

oft X y) - ^-(Bt+y^(x+u)coth((^t/2))+F(x)-F(y))/(2a) 

^ ' ^ ^ 2cJsinh(VIt/2) 



(5-24) X ( Ci{y)I r^-^, ( ^^^^ ^ 

^ ' V ^^^^^/"V(Tsinh(^/At/2)y 

+ C,(,)/„^^/^(— ^g^)). 
As usual F'{x) = f{x)/x and I^diz) = Kd{z) if d is an integer. 

Proof. The proof is similar to the previous cases, with the PDE having 
the same infinitesimal symmetries and invariants. □ 



Again, we will normally have Ci = 1, C2 = 0. 
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Example 5.3. We consider a squared Bessel process of dimensions n. 
We apply the previous result with A = Al? , B = C = Av — 2n+\in? , Ci = I, 

C2 = 0. The PDE ut = 2xuxx + — {^x + has a fundamental solution 

„a ^ = ^ ^-6(x+.)/(2tanh(M)) (y_\'"'^^'\ ^ ( Wxy_\ 

' '^^ 2sinh(6t) \x) V(«-2)2+8z./2l^sinh(6t)y'- 

This reduces to the necessary density and we find 

e^^yp{t,x,y)dy 



r(/3) (cosh 6i + (2A/6) sinh htf 
iP'x cscYiiht) 



26cosh(6t) + 4Asinh(6t) / ' 



where a = sj{n - 2)'^ + 81/, (5 = \{2 + a + n),-f = \{2 + a-n), a = \{a + n + 2) 
and (3=^. 

5.5. An example in the 7 = case. We prove the following result. 



Theorem 5.11. Let f he a solution of the Riccati equation af + ^/^ 

C I B : If A ie,,^\^2 



+ f + ^{A - 16/xcr)x^ A > 0. Then the PDE 

(5.25) Ut = cruxx + f{x)ux - /xx^u, x > 0, 

has a fundamental solution of the form 

\fAxy 



p{t,x,y) 



2smh{VAt/2) 
X exp 



(5.26) 



Bt F{y)-F{x) _ ^A(x2j-y2)^\ 
Aa 2a 8o-tanh(VAt/2) / 

y/Axy \ 



X (Ci{y)Ip 
+ Ci{y)L 

Here F' = f , (3 = + 2C /a'^ and we interpret I-p to he Kp if (3 
integer. 



4cJsinh(VAt/2)/ 

4c7sinh(^/At/2)/ 



IS an 
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Proof. Applying Lie's algorithm shows that if / satisfies the given Ric- 
cati equation, then the finite dimensional part of the Lie algebra of symme- 
tries is spanned by vi = (9^, V4 = udu, 



V2 



8a 4cj 



m = ^~^^/^ ^ fi = ^ 4"^^°^ • Proceeding as previously, we find that the fun- 
damental solution has the form 

v(t X y) = r-'^Bt+VA(x2W)coth(VAt/2)+ialn(smh{VAt/2))/{8a) f ^ \ 

^ ' ^ ^ Vsinh(^/It/2)/ 
Where 16r]^a^v" (rj) - {Ay'^rf + 8C)v{r]) = 0. The result follows. □ 

As previously, we will normally take Ci = 1, C2 = 0. 

Example 5.4. We consider a radial Ornstein-Uhlenbeck process. The 
SDE is dXt = + bXt) + V^dWt, Xq = x. For simpUcity we take a > 1/2. 
The fundamental solution is found by the previous theorem to be 



sinh(at) 

(5.27) 

i.t-a(x2+j;2)/(4tanh(at))-(V4)(x2-?;2) j ( __axy_\ 



a = + 4/Li and = ^(a -|- 1). Note that as /U ^ 0, this reduces to the 
transition density for a radial Ornstein-Uhlenbeck process. (See [9] and [10].) 
From which we find for example that 

£;,(e-^^*"^/o^^''^^) = re-'y'pit,x,y)dy 



2g-fea;2/4+a(a-(6-4A)coth(ia))x2/(4(b-4A-acoth(ta)))-i/t 

(6-4A)sinh(ta)~ 



X ( cosh(ta) 



a 



6. Computable functionals for a given drift. Observe that the roles of 
g{x) and f{x) can be swapped in the sense that we can allow f{x) to be a 
fixed drift and we can then determine the functions g for which the PDEs 
(1.3) has a nontrivial Lie group of symmetries. Thus, even though we cannot 
determine the transition density for all drifts /, we can always compute 
specific functionals (,x.g{Xt) for a given drift. 



30 



M. CRADDOCK AND K. A. LENNOX 



Example 6.1. A simple example illustrates the essential idea. Consider 
the drift /(x) = a — h^/x. Suppose o" = 7 = 1. We take B to be arbitrary, 
but for simplicity we assume that ^ > 0, i? > 0. If xf — / + 1 /2/^ + 2xg{x) = 
Ax + B then clearly the function g must be of the form 

, , A-b'^/2 a-o?/2 + B ab-b/2 

^(^) = ^^ + — i — + ^7^- 

From the stationary solution uo{x) = a;(-^~")/^e^^'/^/ ^ ^^r,^ {'\/2Ax) we find 
that 



-.N rr ^ N x^/^^"/^ ri, X(x + At^ 2)]^ f V2A^ 
(6.1 Ux(x,t) = — exp<^-6x ^ U /nymi 7^ — 



is a solution of the PDE ut = xu^x + (a — hy/x)ux — g{x)u. Inverting the 
Laplace transform to obtain the fundamental solution gives 



p{t,x,y) 

(6.2) 



l/a;x{i-a)/2 ^(x- x + y \ 



X I 



2Jxy 



One can check that this is the correct fundamental solution and so that if 
the diffusion Xt satisfies the SDE dXt = (a - by/Xl) dt + y/2XtdWt, then 

ii;,(e^^^'+/o9(^»)'^^) 

= 2'^(2xt)('^+^-")/2e{-^t+(fe-2/t)^)/2(2 + t{b + 2A))("^^''^'')/2 

r((l + a + ^l)/2) (l + a+jx 2x 
-i-ri ,1 



r(l+/i) ' 'V 2 ' ^'t(2 + i(6 + 2A)) 



where ^ = V1 + 2B. This is not enough to deduce the transition density 
except in the case a = ^. However these calculations can be carried out for 
any one dimensional diffusion and can be potentially quite useful. 
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